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^ ; ABSTRACT 

On. Recent work has shown that at high redshift, the relative velocity between dark matter and baryonic 

gas is typically supersonic. This relative velocity suppresses the formation of the earliest baryonic 
| structures like minihalos, and the suppression is modulated on large scales. This effect imprints a 

(f) . characteristic shape in the clustering power spectrum of the earliest structures, with significant power 

on ^100 Mpc scales featuring highly pronounced baryon acoustic oscillations. The amplitude of these 
oscillations is orders of magnitude larger at z ~ 20 than previously expected. This characteristic 
signature can allow us to distinguish the effects of minihalos on intergalactic gas at times preceding 
and during reionization. We illustrate this effect with the example of 21 cm emission and absorption 
| from redshifts during and before reionization. This effect can potentially allow us to probe physics on 

kpc scales using observations on 100 Mpc scales. 
*-£h ■ We present sensitivity forecasts for FAST and Arecibo. Depending on parameters, this enhanced 

»pHj structure may be detectable by Arecibo at z ~ 15 — 20, and with appropriate instrumentation FAST 

q , could measure the BAO power spectrum with high precision. In principle, this effect could also pose 

a serious challenge for efforts to constrain dark energy using observations of the BAO feature at low 
' redshift. 

: 

1. INTRODUCTION 

Structure formation in the standard inflationary A-Cold Dark Matter (ACDM) cosmological model is expected to 
proceed hierarchically. The earliest bound, virialized structures arise on small scales, and as time progresses these 
■ objects merge and accrete mass, growing ever larger until cosmic acceleration at low redshift freezes out the growth of 
large-scale structure. 

The role played by the earliest generations of collapsed structures in the thermal history of the universe is, at present, 
Q\ • unclear. Observations of the spectra of high-redshift quasars indicate that the reionization of the intergalactic medium 
f^*) ' was largely complete by redshift z w 7 (e.g. Fan et al. 2006), while measurements of the Thomson scattering optical 
, depth of the cosmic microwave background suggest that reionization occurred at z ~ 10 (Larson et al. 2010). At 
,— I ■ these epochs, the typical masses of collapsed dark matter halos range from rare 10 9 M Q objects, down to (plausibly) 
| Earth-mass halos (Green et al. 2005). As we discuss below, the smallest dark matter halos are unable to attract 
•rH . baryons, resulting in an effective lower mass limit near 1O 5 M . 

Even if low-mass halos are able to acquire baryons, they may be unable to convert those baryons into stars (Tegmark 
5h ] et al. 1997; Bromm et al. 2009; Loeb 2010), since star formation requires the presence of cold, dense gas. Objects 
massive enough to attract gas but whose virial temperatures are below ~ 10 4 K, termed minihalos, cannot cool their 
gas through atomic lines and must therefore rely upon molecular cooling processes. It is unclear whether molecular 
processes can cool minihalo gas sufficiently to allow star formation. Prior to the formation of the first stars in the 
Universe, the formation of molecular H2 catalyzed by residual free electrons left over after recombination appears 
insufficient to allow efficient cooling at redshifts z < 20 (Tegmark et al. 1997). However, feedback from the first 
luminous objects can change this result. Positive feedback, for example from ionizing X-rays that strip electrons from 
atoms and thereby spur molecule creation, could lead to efficient cooling. Conversely, negative feedback in the form 
of ultraviolet radiation in the Lyman and Werner bands could destroy molecules and suppress H2 cooling over large 
volumes (Yoshida et al. 2007). Given this wide range of possible scenarios, it is unclear whether minihalos can form stars 
and whether they might be important during the reionization of the intergalactic medium. This uncertainty, however, 
may be viewed as an opportunity: any probe that can quantify the importance of minihalos during redshifts preceding 
and during reionization would dramatically help to elucidate the physics of star formation in the hrst structures that 
arise in the Universe. 

Recently, Tseliakhovich & Hirata (2010, hereafter TH) pointed out an important effect governing the formation of 
~ 1O 5 M minihalos, that had previously been overlooked. As we discuss below, this effect can provide a minihalo 
signature in many potential observables. TH noticed that the relative velocity between dark matter and baryons 
following recombination is typically supersonic. This relative velocity arises because dark matter is accelerated by 



2 



gravitational potential gradients, while baryons are Jeans stabilized against gravitational collapse due to their tight 
coupling with the photon radiation field, until recombination. Because DM and baryons suffer different accelerations, 
they acquire significant relative velocities that are predominantly sourced by potential fluctuations on scales of order 
the sound horizon at recombination, ~ 150 Mpc. At recombination, the baryon sound speed and Jeans length fall 
precipitously, allowing baryons to respond to the same gravitational potential wells that accelerate dark matter. 
Subsequently, the large-scale relative velocity is unsourced and decays as a -1 following recombination. At redshifts 
z w 1000, the relative DM-baryon motion is highly supersonic, with Mach numbers Ai « 5. The gas sound speed does 
not initially decay as quickly as a -1 because of residual thermal coupling to the CMB, so the Mach number diminishes 
over time to M. ~ 2 at z ~ 100, and remains nearly constant thereafter. 

Because the relative motions between baryons and DM are supersonic, they have significant effects on the growth 
of structure. TH computed the effects of relative DM-baryon velocities at high redshift, in the perturbative regime of 
structure formation. They showed that these motions cause a ~ 10% suppression of the matter power spectrum at 
k w 200Mpc _1 compared to standard linear perturbation theory calculations. Using a Press-Schechter approach, TH 
also suggested that the abundance of dark matter halos of mass M ~ 10 6 M Q could be suppressed by a factor of ~ 2 
at redshift z = 40. At later redshifts closer to reionization, z w 10 — 20, the halo abundance would be much closer to 
standard predictions. 

In this paper, we consider a similar effect of the supersonic relative velocities that TH discussed. Instead of considering 
the effect on dark matter halos, we examine the impact on baryonic objects, which likely determine the properties 
of observable quantities like 21 cm absorption, emission, etc. The supersonic flow changes the mass threshold of 
baryonically populated halos, and this effect can be exponentially large for rare objects. The relative bulk flows 
between DM and baryons are modulated on large scales, of order ~ 100 Mpc, meaning that the minihalo abundance 
is similarly modulated. This provides a large-scale signature of the effect of minihalos. 

All calculations presented here assume WMAP (Larson et al. 2010) cosmological parameters: Q TO = 0.27, 17a = 0.73, 
n b = 0.045, h = 0.7, n s = 0.96, a$ = 0.8. 

2. BARYONIC COLLAPSE FRACTION 

In this section, we compute the statistics of the collapsed baryonic density, which we parametrize using f c , the 
collapsed fraction of baryons. The baryonic collapsed fraction is different than the CDM collapsed fraction because 
baryons, unlike CDM, have a nonzero temperature, and therefore cannot fall into shallow potential wells. Naively, we 
might expect that the characteristic mass M c unable to attract baryons would scale like the Jeans mass Mj oc c^/GH, 
however cosmological hydrodynamic simulations (Gnedin 2000; Naoz et al. 2010) have found instead that M c scales 
like the filter scale (Gnedin & Hui 1998), 

M y* = l r da'M 2 J\a') 
a Jo 

For a standard WMAP cosmology, M c ~2-3x 10 4 M Q at z ~ 20 (Naoz et al. 2010). 

Just as gas pressure impedes the collapse of baryons, we expect that bulk relative velocity between gas and dark 
matter halos will also suppress the accretion of baryons. We can make a simple estimate of the minimum halo mass 
that can accrete gas moving at some bulk velocity v relative to the halo through analogy to the above argument. 
Assuming that bulk kinetic energy is converted into thermal energy when gas falls into the halo, we expect that Eqn. 
(1) will still hold if we replace the sound speed c s with c S;e ff = (c 2 + w 2 ) 1 / 2 . Since c s and v both scale as a -1 (see 
below), this effectively multiplies the critical mass scale by the factor (1 + v 2 /c 2 ) 3 / 2 . 

Given M c , we can compute the baryon collapsed fraction f c by integrating the halo mass function, 

/c =^£ M ^ rfM ' (2) 

where we use the fitting function of Sheth et al. (2001) to describe the halo mass function dn/dM. Since M c is a 
function of the local relative velocity between CDM and baryons, v c b = v c — Vb, we see that the collapse fraction f c 
is also modulated by this velocity, f c (x) = f c (v c b(xj). We have assumed that the velocity dependence of the collapse 
fraction is a sharp cut-off in the mass function, whereas hydrodynamic simulations (e.g. Gnedin 2000; Naoz et al. 
2010) find a smooth transition in baryon content of halos below and above M c . This distinction will not be significant 
for our results: the important point is that the collapse fraction is now a function of the local relative bulk velocity 
between dark matter and baryons. 
The relative velocity between CDM and baryons can be computed using the linearized continuity equation 

<5 + a _1 V • v = 0, (3) 

using comoving coordinates rather than proper coordinates. In the linear regime, where we have potential flow, the 
velocity is 
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Fig. 1. — (Left) The power spectrum of v c ^ fluctuations, A^(fc) = P pl i(k)[aT(k)] 2 /2n 2 k, at z = 15. (Right) Velocity correlation functions 
■01 (blue), ip r = ipi -+- ip2 (red), and ip v = ip\ + ip 2 /2 (black). Dotted curves depict regions where the correlation function is negative. 
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For adiabatic perturbations, CDM and baryons have the same primordial overdensity perturbations S pri , so the relative 
velocity is 

.ak , 



v cb (k,a) = -i—T cb (k,a)S pli (k), 



(5) 



where T cb — T c ~ T b is the difference between the CDM and baryon transfer functions. The two-point correlation 
function for velocities is then 
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Here, P pr i(k) is the primordial density power spectrum, (#pri(fei)#pri(fe2)) = (27r) 3 P pr i(fci)^' 3 ' (fei + £2). Note that 
fpri(fc) = ^4 fc™ s is time-independent; all of the time dependence of the power spectrum is contained in the transfer 
function. We calculate the time-dependent CDM and baryon transfer functions using CAMB (Lewis et al. 2000). Figure 
1 illustrates these velocity correlations. On small scales, velocities at nearby points are almost perfectly correlated, 
ipl ~ 1 and ip2 ~ 0. Towards larger scales of order the sound horizon ~ 150 Mpc, the radial component of the velocity 
is anti-correlated, since the relative CDM-baryon velocities are sourced by structures on these scales. The correlations 
fall off steeply on scales much larger than the sound horizon. 

Since we have assumed that the collapse fraction f c is a function of the local (relative) velocity v, we can easily 
compute the probability distribution of f c and its moments by integrating over the Gaussian distribution of v. For 
example, the mean collapse fraction is 
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P(v)f c (v)d 3 v 
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We plot the redshift dependence of (/ c ) in Fig. 2. 
Similarly, the two-point correlation function of / c is 



(f c (x)f c {x + r)) - / d 3 v 1 d 3 v 2 P(v 1 ,v 2 )f c (v 1 )f c (v 2 ) 
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Fig. 2. — Redshift dependence of the mean collapse fraction (/ c ) (solid black curve), the collapse fraction in the absence of relative 
baryon-DM motions f c (v = 0) (dotted black curve), and the square of the effective bias, (dashed blue curve), defined in Eqn. (14). 

where v\ and v 2 are the velocities at points x and x + r, respectively, whose probability P is a Gaussian distribution 
with covariance given by Eqn. (6). This is a 6-D integral, however two of the integrals are elementary. If we write 
u = v I <T\ with components u r and Ut parallel and perpendicular to the separation vector, then 



(f c (x)f c (x + r)) = 
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where ipr — ipi + "02 and Iq(x) is a modified Bessel function. In various regimes, this integral may be simplified by 
Taylor expansion. For example, at large radius where |^i| <C 1 and |^v| <C 1, the two-point function is approximately 



(f c (x)f c (x + r))^- / du 1 du*u\ulf c {v 1 )Uv 2 )e-^+<)! 2 
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Writing 1 + C/(r) = (f c (x)f c (x + r))/(f c ) 2 , we find 
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(14) 



where a 1 = (v 2 ) = 3cr 2 is the 3-D velocity dispersion. 

This expression has a straightforward interpretation. We have assumed that the collapse fraction is a function of the 
local velocity v. Therefore, on large scales, f c will be a biased tracer of v 2 . The two-point correlation function of v 2 is 
simply 2ip 2 + if; 2 . So Eqn. (14) is not surprising: f c is indeed a biased tracer of v 2 , with a large-scale bias coefficient 
frctr = ("f 2 /c)/c 2 (/c) — 1. The form of this result is therefore generic. The numerical value of the bias coefficient 
depends on the precise details of the baryon content of low-mass halos, which we have assumed to be a step-function 
at M c for simplicity. A more realistic model will have a slightly different bias coefficient, but the result that £/ oc ip v 
holds generically as long as the gas content of halos depends on the local relative bulk velocity. 

In Figure 3 we plot the two-point correlation function and power spectrum of the baryonic collapse fraction. The 
two-point function exhibits strong clustering on large scales, roughly ~ 10% fractional fluctuations on 100 Mpc scales 
at z = 20. There is also a pronounced baryon acoustic oscillation feature, which is more plainly visible in the power 
spectrum. For comparison, we have also plotted the matter power spectrum at the same redshift. Two aspects to 
note are that the collapse fraction has considerably enhanced power on large scales, compared to the matter power 
spectrum, and that the amplitude of the baryon oscillations are far larger in the f c power spectrum. As the figure 
illustrates, the baryon wiggles are typically percent level in the matter P(k), while in the f c power spectrum, the 
baryon oscillations are order unity. 
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Fig. 3. — (Left) The f c correlation function §/(»") at z = 10 (red), z = 20 (blue) and z = 30 (black). The solid curves show the full 
correlation function calculated from Eqn. (12), while the dotted curves show the linear bias approximation of Eqn. (14). At small radii, 
the linear bias approximation begins to break down, when Poisson-like clustering can become significant. (Right) The dimensionless power 
spectrum A 2 = fc 3 P(fc)/27r 2 for fractional fluctuations in collapse fraction 8f = fc/(fc) — 1 (black) and matter <5 = p/p — 1 (blue) at z = 20. 

The large-scale modulations of collapsed baryonic mass found in the previous section can have several observable 
consequences. In this section, we discuss one such consequence: the effect of relative baryon-DM velocities on the 21-cm 
anisotropies preceding reionization. The physics of 21-cm fluctuations at high rcdshift has been discussed extensively 
in the literature; see Madau et al. (1997) and Furlanetto et al. (2006) for reviews. 

Prior to reionization, the gas kinetic temperature X^m is far below the CMB temperature: roughly speaking, T^ n ~ 
T cia b x (1 + z)/150 for z < 100. When the first stars begin to emit Lyman-a radiation, the Wouthuysen-Field effect 
couples the spin temperature of neutral hydrogen atoms to the gas kinetic temperature, allowing neutral gas to be seen 
in 21-cm absorption against the cosmic microwave background. Neglecting atomic collisions, the spin temperature is 
given by 
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where the coupling coefficient x a may be expressed in terms of the ratio of Lyman a photons to H atoms, n a , as 
(Furlanetto et al. 2006) 
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Given the spin temperature T s , the optical depth to 21 cm absorption is (Madau et al. 1997) 
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and the observed (i.e. redshifted) 21 cm brightness temperature contrast against the CMB is then 
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(16) 
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where Tqmb(z) = 2.726(1 + z) K. Therefore, given the Lyman-a intensity at any given point, we can compute the 
observed temperature contrast. 

We estimate the Lyman a intensity by assuming that, on average, each collapsed baryon emits N a Lyman a photons 
that escape into the intergalactic medium. For simplicity, we assume that photons are emitted with a flat spectrum, 
vdN I dv = const, but it is straightforward to generalize our results for an arbitrary emission spectrum. The local 
Lyman-a photon number density at each point in space and redshift is then given by a convolution of the collapsed 
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baryon density with a retarded Green's function, 

„._ w „/A*$(,,,)*!a£M, (21) 

or in Fourier space, 

f'° df 

n a (k,r] Q )=N a dr]^(k,r])j Q (kc(r]o-v)) (22) 
Jo dr) 

where rj is conformal time and jo(x) — sm(x)/x. We can further simplify this expression by using our previous result 
that, on large scales, f c (k) has time dependence oc (f c )b c e = (fc) — (v 2 f c )/(J 2 . Therefore, we have 

n a (k,z) =N a f c (k,z)W(k,z) (23) 

where the smoothing filter is 

W ^ z ^ = tt^t r d v^r^MMvo »?)). (24) 

(jc)Ocff Jo dr] 

In this expression, 770 is the conformal time at redshift z, and note that d/drj — —Hd/dz. 

The propagation of Lyman a photons over large distances considerably damps the spatial fluctuations in Ly a 
intensity. However, our discussion so far has neglected an important effect: Lyman a photons can only travel a limited 
distance. A rest frame Lyman a photon which participates in pumping was emitted at some distance bluewards of 
Lyman a, redshifting as it travels. The higher the frequency at emission, the longer the distance that the photon 
travels before redshifting into Lyman a. However, a photon that is emitted at a wavelength shorter than Lyman 
P will be absorbed in the neutral intergalactic medium, and ultimately lost to double photon decay, before it can 
redshift into Lyman a. This means that gas clouds at redshift z cannot be pumped by photons emitted by sources at 
redshift z C mit > Zhor, where (1 + Zhor) = (32/27) x (1 + z). This gives a natural maximal horizon distance for Lyman a 
pumping. Conceivably, the propagation distance could be even shorter given sufficient molecular opacity, either in the 
host minihalos or in the intergalactic medium (Ricotti et al. 2001), but we disregard this possibility in our calculations. 

If Zh or /z — 1 -c 1, we can approximate Eqn. (24) as 

w{k ' z) ~-m*rJ z dz *^- 3 °w z *- z) )- (25) 

This shows how the shape of the smoothing window depends on the formation history of collapsed baryonic objects 
and their Lyman a emission. As noted above, the shape of W(k, z) also depends on the spectrum of escaping UV 
emission from the first stars, which will generally be much more complicated than we have assumed here. Fortunately, 
it is entirely straightforward to compute how W changes when realistic spectra and opacity are used instead of the 
flat spectrum and sharp cutoff that we have adopted for simplicity. Given this expression for the smoothing window, 
we can compute the number of Lyman a photons per atom, n a , using Eqn. (23), which then gives x a and the spin 
temperature T s using Eqns. (15-16). Given T s , we compute the optical depth and brightness temperature using Eqns. 
(19-20). 

The Lyman-a intensity n a is a linear function of the collapse fraction, so its power spectrum is simply the product of 
the f c power spectrum with the (square of the) window function, Eqn. (24). The brightness temperature is a nonlinear 
but local function of n a . Therefore, on large scales it is a biased tracer of the intensity field, and its power spectrum 
will be proportional to the n a power spectrum, with some proportionality coefficient. We could write down an analytic 
expression for this bias coefficient in terms of the ./V-point correlation functions of n a , but it is simpler to calculate it by 
simulation instead. Accordingly, we have generated realizations of the brightness temperature field. We first generate 
realizations of the Gaussian random relative velocity field v c b, which we then transform into collapse fraction f c using 
Eqns. (1-2), replacing c s — > c SjC ff as described in §2. From the collapse fraction, we compute the Lyman-a intensity 
using Eqn. (23), which then gives the spin temperature T s , optical depth r and brightness temperature contrast 6Tb 
using Eqns. (15-20). We generate realizations in a 2 h~ 1 Gpc box of 1024 3 pixels at a variety of different rcdshifts, 
varying the number of Lyman-a photons per collapsed baryon, N a . 

Figure 4 illustrates the brightness temperature power spectrum. As expected, the 6Tb power spectrum is proportional 
to the product of the f c power spectrum, multiplied by the square of the window function W(k, z) given by Eqn. (24), 
which suppresses small-scale fluctuations in the brightness temperature. The shape of the power spectrum is therefore 
simple to calculate. The power spectrum peaks on the scale of the Lyman-a horizon, and exhibits damped but 
pronounced acoustic oscillations at higher wavevectors. 

The amplitude of the power spectrum is a nontrivial function of N a and z. At high redshift, when the collapse fraction 
is small and the Lyman-a pumping intensity n a is weak, the spin temperature is close to the CMB temperature, with 
small fluctuations proportional to n a . The brightness temperature contrast therefore grows rapidly in time. Eventually, 
however, the spin temperature begins to saturate at the gas kinetic temperature 7ki n - As n a becomes very large, the 
spin temperature begins to approach a uniform value everywhere, T s — ► Tk; n , and so the brightness temperature 
fluctuations actually diminish with increasing n a . Accordingly, the 6Tb power spectrum peaks and then decreases 
towards lower redshift. Of course, our assumption that T^ n remains close to its adiabatic value becomes suspect in 
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Fig. 4. — (Left) The 21 cm temperature power spectrum at redshift z = 20, for N a = 100 Lyman a photons per collapsed atom. The 
blue curve with error bars shows the brightness power spectrum from UV sources that trace minihalos, while the red dashed curve shows 
the fluctuations expected for uniform Lyman a radiation that pumps gas with density fluctuations tracing the matter density fluctuations 
at this redshift. Since density and velocity are uncorrelated, these two power spectra add linearly (solid black curve). (Right) RMS 
temperature fluctuations 8T b as a function of redshift z for a Gaussian beam with a diffraction limit of a 500m dish (such as FAST). A 
300m dish measures a slightly lower signal. Different curves are: N a = 10 (red dot-dash), 30 (blue dashed) and 100 (black solid). The 
shaded area indicates the 1 — a errors per unit redshift expected for 1000 days of observations using a 500m dish like FAST. 

this regime. Eventually, a sufficiently strong Lyman-a intensity will not merely pump the 21-cm transition, but will 
appreciably heat the gas as well. 

4. OBSERVATIONAL PROSPECTS 

In order to measure the very large scale structure, one needs high brightness sensitivity on BAO scales. At z > 10, 
the angular scale of BAO changes little, and ~ 10' scales are important, just as they are in the CMB. 

To achieve high brightness sensitivity, a filled aperture is desirable. Telescopes such as Arecibo and FAST (Smits 
et al. 2009) would be well suited. With filled apertures of 300m and 500m respectively, their angular resolution for 21 
cm at z ~ 20 is 40 and 30 arc minutes, respectively. We focus our attention on FAST, which is expected to observe 
at sufficiently long wavelengths. At these low frequencies, sensitivities are sky limited, with T sys ~ 3000K. It is also 
straightforward to observe with a focal plane array, so we reference our forecasts to a 100 pixel array, which would be 
a 40m dipole array. Such a focal plane array would allow primary beam illumination to compensate with frequency to 
make frequency independent beams, enabling accurate foreground subtraction. The maximum transverse k± = 0.08. 

We further assume the sky is drift scanned, which minimizes systematic errors, and makes this experiment comparable 
to the recent GBT 21cm intensity detection (Chang et al. 2010). For a square array, the field of view is 10 beam width, 
or about 5 degrees. We use half the scanned area as useful, allowing for galaxy and point source cuts. 

At zenith for a latitude of 30 degrees, this scans 1000 square degrees. We use an integration time of 1000 days on 
the sky. The exposure time per pixel is 5 x 10 6 seconds. In a bandwidth of 2 MHz, this results in a pixel noise of 1 
mK, well matched to the expected signal. This map contains ~ 10 pixels, so one expects to measure the power on 
the beam scale with ~ lOOer. 

The specific forecast parameters used for the noise estimate in figure 4 are: T s k y = 300(z//150MHz) -2 - 7 , 24 beams 
on the sky (the equivalent of a LOFAR station signal processor), and a diffraction limited beam. 

5. DISCUSSION 

We have investigated the effect of relative motions between baryons and dark matter on the formation of the smallest 
galaxies at high redshift, z ~ 20. The formation of galaxies in minihalos of mass M < 10 6 M Q is modulated by large- 
scale bulk velocities between gas and dark matter, and so the clustering of these objects will contain a contribution 
proportional to the relative velocity two-point function. The velocity power spectrum exhibits significant correlations 
on large scales of order 100 Mpc, with pronounced baryon acoustic oscillations. Accordingly, the large-scale clustering 
of minihalos exhibits similar features, and the large-scale correlations of any observable that traces minihalos will show 
similar behavior. We illustrated this in the previous section with a calculation of the 21 cm absorption power spectrum 
prior to reionization, when the kinetic temperature of intergalactic gas is much colder than the CMB temperature. 
Our predictions for 21 cm absorption are model dependent, principally depending on the number of Lyman a photons 
emitted per collapsed baryon. 

A similar argument holds for any observable that traces minihalos. For example, suppose that minihalos make 
a significant contribution to the ionizing flux at the time of reionization. Then the large-scale power spectrum of 
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Fig. 5. — The 21 cm temperature power spectrum AT 2 = fc 3 P(fc)/2-7r 2 at z = 10 (solid black), assuming that the ionized fraction is 
proportional to the local baryon collapsed fraction, x e oc f c . For comparison, we also plot the matter power spectrum multiplied by 
(1 + z) X (9.4mK) 2 , as the dotted blue curve. 

patchy reionization will contain a contribution proportional to the bulk velocity power spectrum, as can easily be seen. 
During reionization, let us assume that T s 3> Tomb, and that fluctuations in the ionized fraction x e trace the local 
collapsed fraction f c . In this regime, the observed temperature contrast in Eqn. (20) becomes independent of the spin 
temperature 

ST b ~ 0.155K ft,— 7^(1 + z) 1/2 x e w 9.4mK (1 + z) 1/2 x e . (26) 

We fix the proportionality constant relating x e and f c by requiring that the Thomson scattering optical depth r cannot 
exceed the WMAP bound, r = 0.088. We assume that reionization completes at z IC — 7, and assume that x e at z > z Te 
is patchy and traces f c . This provides an upper bound on the patchy reionization signal from minihalos. The resulting 
power spectrum is shown in Fig. 5. As expected, the minihalo contribution to the power spectrum shows pronounced 
baryon acoustic oscillations, which will be at least partially smoothed when we account for the ~ 5 — 10 Mpc sizes of 
ionized bubbles during reionization. 

Although this signal could be present, detecting it during reionization may be challenging in practice, since there are 
other sources of large-scale power in 21-cm correlations. For example, the clustering of the ionizing sources will reflect 
the clustering of the halos hosting those sources, which trace the clustering of the matter field. Because matter density 
fluctuations are uncorrelated with bulk velocity fluctuations, these two terms simply add in the overall power spectrum 
of the sources. Realistically, we expect the clustering from bulk velocities to be subdominant to the clustering from 
density perturbations over many of the scales of interest during reionization (see Fig. 5). On the other hand, TH 
have shown that the halo power spectrum itself also contains a term proportional to the relative bulk velocity power 
spectrum, for halos of order M ~ 10 6 M Q . 

Another place where these effects may be observable is in the CMB itself. At small angular scales, the so called 
kinetic Sunyaev-Zeldovich effect, caused by radial Doppler motions of electrons off which photons scatter, can be 
important (e.g. Hernandez-Monteagudo & Ho 2009). This effect is enhanced if reionization is patchy (McQuinn et al. 
2005). If these reionization patches have suppressed bulk velocities then the effect will be reduced. On the other 
hand, the large scale correlation of reionization patches predicted here will lead to large scale fluctuations in optical 
depth, inducing a modulated suppression of primary CMB fluctuations, which may be detectable with a higher order 
correlation analysis in CMB. 

At low redshifts, z < 5, the typical masses of collapsing halos are generally much larger than minihalo masses, and so 
we would naively expect the effects of relative velocities between DM and baryons to become unimportant. However, 
in principle the signatures of minihalos could persist even in late-time observables. For example, if reionization is 
patchy on ~ 100 Mpc scales due to minihalos, the subsequent star formation history inside patches that reionize early 
could differ from patches that reionize later. This could lead to spatial variations in galaxy formation at late times, 
on scales of order the BAO scale. This is potentially worrisome for BAO probes of dark energy (e.g. Schlegel et al. 
2009), which rely on precise determination of the BAO feature in the galaxy two-point function. Any contamination 
from minihalo effects could shift the location of the BAO feature and thereby create a bias in measurements of the 
equation of state parameter w, analogous to the results of Pritchard et al. (2007). Removing this source of systematic 
uncertainty appears daunting. Eisenstein et al. (2007) have shown that BAO probes of dark energy are quite insensitive 
to smooth distortions to the shape of the power spectrum. However, the minihalo effect could be far more pernicious, 
since the velocity power spectrum is not smooth, but has pronounced baryon oscillations that are presumably out of 
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phase with the oscillations in the matter power spectrum. Marginalizing over an unknown minihalo contaminant could 
significantly degrade BAO constraints on dark energy. This underscores the need for more work on this subject. 

In summary, the supersonic relative motions between baryons and dark matter can dramatically affect the formation 
of the earliest collapsed baryonic objects in the Universe. We expect the abundance of minihalos of mass M ~ 10 5 M© 
to be modulated on ~ 100 Mpc scales, and any tracer of minihalos should show the same modulations, with highly 
pronounced baryon acoustic oscillations. Specific predictions for any observable are necessarily model dependent, but 
for plausible scenarios this large scale signal is within the sensitivity range of existing and upcoming observatories. 
This effect could allow us to detect the signature of the earliest galaxies, written across the sky. 

We thank Chris Hirata for useful discussions. The calculations presented in this paper have made use of publicly 
available software, including CAMB 1 , CUBA 2 , FFTLog 3 and GSL 4 . We thank the authors of these libraries for making 
their software public. This work is supported by the Swiss National Foundation under contract 200021-116696/1 and 
WCU grant R32-2009-000-10130-0. 
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